Astronomy &: Astrophysics manuscript no. chauvin'dynamicsHdl96885 


© ESO 2010 


September 30, 2010 





Planetary systems in close binary stars: the case of HD 196885 

Combined astrometric and radial velocity study 

G. Chauvin 1 , H. Beust 1 , A.-M. Lagrange 1 , and A. Eggenberger 1 
^aboratoire d'Astrophysique, Observatoire de Grenoble, UJF, CNRS; BP 53, F-38041 GRENOBLE Cedex 9 (France) 
Received ; accepted 

ABSTRACT 

Context. More than fifty candidate planets are presently known to orbit one component of a binary or a multiple star 
system. Planets can therefore form and survive in such an environment, although recent observing surveys indicate that 
short-separation binaries do not favour the presence of a planetary system around one of the component. Dynamical 
interactions with the secondary component can actually significantly impact the giant planet formation and evolution. 
For this reason, rare close binaries hosting giant planets offer an ideal laboratory to explore the properties and the 
stability of such extreme planetary systems. 

Aims. In the course of our CFHT and VLT coronographic imaging survey dedicated to the search for faint companions of 
exoplanet host stars, a close (~ 20 AU) secondary stellar companion to the exoplanet host HD196885 A was discovered. 
In this study, our aim is to monitor the orbital motion of the binary companion. Combining radial velocity and high 
contrast imaging observations, we aim to derive the orbital properties of the complete system and to test its dynamical 
stability to reveal its formation. 

Methods. For more than 4 years, we have used the NaCo near-infrared adaptive optics instrument to monitor the 
astrometric position of HD 196885 B relative to A. The system was observed at five different epochs from August 2005 
to August 2009 and accurate relative positions were determined. 

Results. Our observations fully reject the stationary background hypothesis for HD196885B. The two components are 
found to be comoving. The orbital motion of HD 196885 B is well resolved and the orbital curvature is even detected. From 
our imaging data combined with published radial velocity measurements, we refine the complete orbital parameters of 
the stellar component. We derive for the first time its orbital inclination and its accurate mass. We find also solutions for 
the inner giant planet HD196885 Ab compatible with previous independent radial velocity studies. Finally, we investigate 
the stability of the inner giant planet HD196885 Ab due to the binary companion proximity. Our dynamical simulations 
show that the system is currently and surprisingly more stable in a high mutual inclination configuration that falls in 
the Kozai resonance regime. If confirmed, this system would constitute one of the most compact non-coplanar systems 
known so far. It would raise several questions about its formation and stability. 
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1. Introduction 

Among current exoplanets hunting techniques, radial ve- 
locity (RV) measurements are nowadays the most success- 
ful method for detecting exo-planetary systems (Udry & 
Santos 2007; Cumming et al. 2008). Although the min- 
imum mass is derived, a small correction including 
the distribution of inclination is expected to access 
the true mass distribution of the current candidate 
planets (Watson et al. 2010; Jorrissen et al. 2001). 
Originally focused on quiet solar-type stars that show nu- 
merous thin absorption lines, RV surveys have recently di- 
versified their samples to consider a broader class of pri- 
mary stars. Telluric planets arc now preferentially searched 
for around M dwarfs as the habitable zone can be ex- 
plored (Mayor et al. 2009; Charbonneau et al. 2009). Giant 
planets have probably been discovered around intermediate 
mass objects (Lagrange et al. 2009), giant stars (Dollingcr 
et al. 2009; Lovis et al. 2007) and have actively been for 
earched around young active stars (Setiawan et al. 2008). 

Send offprint requests to: G. Chauvin 



Finally, planets are now scrutinized in multiple stellar sys- 
tems where RV surveys used to exclude them. This was 
mainly related to the difficulty of processing the multiple 
stellar RV signals at the required precision to detect plan- 
ets (Konacki 2005; Eggenberger & Udry 2007; Toyota et al. 
2009; Desidera et al. 2010). Binaries and triple are particu- 
larly interesting to test the predictions of the planetary for- 
mation and evolution processes. They enable us to under- 
stand how a perturber will impact the planetary system for- 
mation and dynamical evolution. Originally, planets were 
rapidly found in close spectroscopic binaries (Gliese 86, 
Qucloz et al. 2000; 7 Cep Hatzes et al. 2003), confirm- 
ing that circumstellar planetary systems could form and 
survive in such a hostile environment (binary separation 
< 20 AU). Circumbinary planets, undiscovered up until 
now, remain more difficult to detect due to the limited size 
of the sample and only the recent start-up of dedicated 
programmes. 

A few years ago, studies and surveys had success in tack- 
ling the problem of duplicity in planetary systems. Zucker 
& Mazeh (2002) initially showed that the planet properties 
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in binaries were different to planets around single stars. 
Eggenberger et al. (2004) confirmed that the most mas- 
sive short period planets were found in binary systems and 
that planets with short (P < 40 days) periods in binaries 
were likely to have low eccentricities. Tamuz et al. (2008) 
finally found that extremely high eccentric planets were all 
in binaries. In parallel to Doppler surveys, imaging sur- 
veys studied the effect of duplicity on the planet occurence 
around longer period (> 100 AU) binaries. Catalogs com- 
piled by Raghavan et al. (2006), Desidera & Barbieri (2007) 
and Bonavita & Desidera (2007) did not find any signif- 
icant discrepancies in the planet frequency between bi- 
nary and single stellar systems. Only high contrast imag- 
ing studies have been able to access intermediate-separation 
(20 — 100 AU) binaries where the companion influence is the 
most expected. Several deep imaging studies have revealed 
a number of additional close companions to known exo- 
planetary hosts (Patience et al. 2002; Chauvin et al. 2006; 
Eggenberger et al. 2007; Mugrauer et al. 2009). However, a 
dedicated survey using a reference sample was mandatory 
to conduct a proper multiplicity study to test the impact 
of duplicity on the giant planet occurence. Eggenberger et 
al. (2007) observed two subsamples of more than 50 stars 
each. They found a lower corrected binary fraction for the 
planet-host subsample, particularly for physical separations 
shorter than 100 AU (Eggenberger et al. 2008). This sta- 
tistical result corroborates the theoretical prediction that 
multiplicity has a negative impact on planetary formation 
and evolution at less than 100 AU (e.g. Nelson et al. 2000; 
Mayer et al. 2005; Thebault et al. 2006). 

Binaries with small semi-major axes (< 20 AU) 
hosting a planetary system offer an ideal opportu- 
nity to characterize in more detail the dynamical impact 
of the binary companion on the inner circumstellar planet. 
They can also serve as a testbed for planet formation mod- 
els, as they push many of these models parameters to their 
limits. We conducted such a study combining RV and AO 
imaging observations for the Gliese 86 system (Lagrange et 
al. 2006). We confirmed that the companion was a white 
dwarf identified in spectrocopy by Neuhauser et al. (2005). 
We showed that the planetary system around Gliese 86 A 
had survived the later stages of evolutiion of the white 
dwarf progenitor, a probable late-F to early-K type star (i.e 
the mass loss of the B component and a semi-major axis 
intially reduced to 13 AU compared with its 18 AU current 
value). In the course of the VLT/CFHT deep imaging sur- 
veys of exoplanet hosts (Chauvin et al. 2006, 2007), we dis- 
covered a close (~ 20 AU) binary companion to HD 196885 
(see Fig. 1). This system is well suited for a detailed dynam- 
ical study. In this paper, we summarize the observations 
and the results of our recent imaging campaign aimed at 
monitoring the binary orbital motion. We report the best 
orbital adjustment to fit the combined imaging and RV ob- 
servations of the complete system. Finally, we discuss the 
dynamical stability of the inner planetary system that has 
formed and survived despite the dynamical influence of the 
close binary companion. 

2. The HD196885 exoplanet host binary 

HD196885 is an F8V (V = 6.398, B -V = 0.559) star 
located at 33.0 ± 0.9 pc (Perryman et al. 1997). Based 
on CORALIE spectra, Sousa et al. (2006) derived a spec- 
trocopic temperature, surface gravity and metallicity of 
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Fig. 1. VLT/NACO image obtained in i^ s -band using the 
S27 CONICA platescale. The binary is clearly resolved. 

T off = 6340 ± 39 K, log(g) = 4.46 ± 0.02 and [Fe/H] = 
0.29 ± 0.05 respectively. The star v sini was estimated to 
7.3±1.5 kms" 1 from ELODIE spectra (Correia et al. 2008). 
These results were recently supported by independent Lick 
observations reported by Fischer et al. (2009). The chro- 
mospheric activity level is relatively low. Bolometric lumi- 
nosity correction and evolutionary model predictions lead 
to an estimate of the luminosity and the mass of 2.4 L@ 
and 1.3 Mq respectively. Finally, the corresponding stellar 
age derived from evolutionary tracks and from the activ- 
ity level varies between 1.5 to 3.5 Gyr (see Correia et al. 
2008; Fischer et al. 2009). As noted by Correia et al. (2008), 
HD196885 might be part of a wider binary system with 
the star BD+104351 P^] located 192 " north (that would 
correspond to a minimum physical separation of 6330 AU, 
Dommanget et al. 2002). 

In 2004, a RV variation was measured and adjusted with 
a preliminary orbital period of P = 0.95 yrs. The result 
was temporarily reported on the California Planet Search 
Exoplanet Web site and was thus withdrawn due to a sig- 
nificant residual drift in the orbital solution. Nevertheless, 
this star was included in our deep imaging survey of stars 
hosting planets detected by RV observations (Chauvin et 
al. 2006). Our imaging and spectroscopic observations led 
to the detection of a co-moving M1±1V dwarf companion 
located at only 0.7" (23 AU in projected physical separa- 
tion) and likely to be responsible for the trend seen in the 
Lick RV residuals (Chauvin et al. 2007). Our two epochs of 
observations resolved the binary orbital motion confirming 
its physical nature. 

Using a double-Keplerian model for the binary star 
and the planet to adjust their ELODIE, CORALIE and 
CORAVEL observations spread over 14 years, Correia 

1 and not BD+104251B as mentionned by Correia et al. 
(2008) 
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Table 1. Relative astrometric position of the binary companion HD196885B to the cxoplanet host star HD 196885 A. 



UT Date 


Filter/Obj 


Aa 


AS 


AK S 


Plate Scale 


True North 






(mas) 


(mas) 


(mag) 


(mas) 


(o) 


01/08/2005 


K s /B27 


658.5 ± 1.7 


273.6 ± 1.6 


3.08 ±0.05 


27.01 ±0.05 


0.09 ±0.12 


26/08/2006 


K a /$27 


649.8 ± 1.9 


292.6 ± 2.4 


3.06 ±0.04 


27.01 ±0.05 


0.02 ±0.20 


25/08/2007 


KJS27 


640.2 ±3.1 


309.7 ±3.2 


3.07 ±0.06 


27.01 ±0.05 


0.05 ±0.15 


28/06/2008 


K s /S27 


630.9 ±3.1 


326.3 ± 3.2 


3.00 ±0.10 


27.01 ±0.05 


-0.25 ±0.14 


27/08/2009 


KJS27 


614.5 ±3.0 


342.1 ± 2.9 


3.14 ±0.09 


27.01 ±0.05 


-0.35 ±0.08 



et al. (2008) derived a first range of orbital solutions. 
They revised the planet solution with a minimum mass of 
M Ab smi = 2.96 M Jup , a period of P = 3.69 ± 0.03 yrs 
and an eccentricity of e = 0.462 ± 0.026. Moreover, they 
found additional constraints for the binary companion 
HD196885B with a period of P > 40 yr, a semi-major axis 
a > 14 AU and a minimum mass of Mb sini > 0.28 M . 
Based on Lick observations, Fischer et al. (2009) recently 
derived consistent results for both the inner planet and the 
binary companion. Nevertheless, in both studies a relatively 
large range of masses and periods remains for the binary 
companion. Additional observing constraints are therefore 
needed to refine the binary companion properties and to 
understand how it could have affected the formation and 
the stability of the circumstellar planetary system around 
HD 196885 A. 



3. Observation and Data Reduction 

The orbit of HD 196885 AB was monitored with the NACO 
(NAOS-CONICA) high contrast Adaptive Optics (AO) im- 
ager of the VLT-UT4. The NAOS AO system (Rousset et 
al. 2002) is equipped with a tip-tilt mirror, a 185 piezo actu- 
ator deformable mirror and two wavefront sensors (Visible 
and IR). Attached to NAOS, CONIC A (Lenzen et al. 
1998) is the near infrared (1 — 5 /im domain) imaging, 
Lyot coronagraphic, spectroscopic and polarimetric cam- 
era, equipped with a 1024 x 1024 pixels Aladdin InSb ar- 
ray. Observations were obtained at five different epochs in 
August 1st 2005, August 26th 2006, August 25th 2007, 
June 28th 2008 and August 27th 2009. Over the differ- 
ent observing campaigns, the atmospheric conditions were 
sufficiently stable to close the AO loop and resolve both 
components. The offset position of HD 196885 B relative to 
A was well monitored at each epoch. The typical observ- 
ing sequence included a set of five jittered images obtained 
using the K s filter and the S27 camera CONICA (mean 
plate scale of 27.01 mas/pixel). Thsi led to a total expo- 
sure time of ~ 5 min on source. To calibrate the plate scale 
and the detector orientation, we observed the astrometric 
field of flOrilC (McCaughrean & Stauffer 1994). When 
not observable in June 2008, we used as secondary calibra- 
tor the astrometric binary IDS 21506S5133 (van Dessel & 
Sinachopoulos 1993), recalibrated with the 9 Oril C field. 

After cosmetic reductions (dead and hot pixels, dark 
and flat) using eclipse (Devillard 1997), we applied the de- 
convolution algorithm of Veran & Rigaut (1998) to obtain 
the offset position of HD 196885 B relative to A at each 
epoch. Single stars of similar brightness observed on the 
same night were used for point spread function (PSF) esti- 
mation. The results are reported in Table [l] The platescale 
and the true North orientation of the detector are given 
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Fig. 2. VLT/NACO measurements with uncertainties in 
the offset positions of HD 196885 B relative to A, obtained 
in August 1st 2005, August 28th 2006, August 27th 2007, 
June 26th 2008 and August 27th 2009. The solid line 
gives the expected variation of the offset positions 
of B relative to A if B is a background station- 
ary object. This variation takes into account the 
initial offset position of B relative to a in August 
1st 2005 and the proper and parallactic motions of 
HD 196885 A. After April 2006, the predicted offset 
positions go beyond the (Aa, AS) astrometric range 
considered for this figure. 

for each measurement. The relative positions are shown in 
Fig. [2] The expected variation of offset positions, if B is a 
background stationary object, is shown (solid line). It takes 
into account the initial offset position of B relative A and 
the proper ((pL a ,^s) = (47.5 ± 0.9, 83.1 ± 0.5) mas/yr) and 
parallactic (rr = 30.31 ±0.81 mas) motions of HD 196885 A. 
The orbital curvature is detected and the background sta- 
tionary hypothesis for HD 196885 B is fully excluded. 

4. Orbital solution 

In attempt to constrain the physical and orbital proper- 
ties of the HD 196885 system, we considered simultane- 
ously the astrometric data points listed in Table [l] and 
all data from past RV surveys: CORAVEL (9 RV mea- 
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Fig. 3. NaCo astrometric observations and orbital solution of HD 196885 A and B. The full orbital solution is sketched in 
red and superimposed to the astrometric data points in blue. Predicted positions from the fit are shown in green. Left, 
Full orbital solution. Right, Zoomed on with the astometric data points and their uncertainties. 



surements from June 1982 to August 1997), ELODIE (69 
measurements from June 1997 to August 2006), CORALIE 
(33 measurements from April 1999 to November 2002) and 
Lick (75 measurements from 1998 to 2008) spanning over 
26 years (see Correia et al. 2008, Fischer et al. 2009). We 
used the iterative Levenberg-Marquardt x 2 minimization 
method (Press et al. 1992) to simultaneously fit both set 
of data and derive the characteristics of both orbits (the 
outer AB orbit between HD 196885 B and A, and the Ab 
orbit between the giant planet and HD 196885 A). The basic 
assumption is that both orbits contribute to the RV signal 
and that only the AB orbit is seen in the astrometric mo- 
tion. For the RV data, separate zero point velocities have 
been fitted for each set of data to account for instrumental 
shifts. The fitting routine turned out to rapidly converge 
towards a unique best solution reported in Table [2j This 
solution is also sketched in Figs. [3] and [4] and superimposed 
to the astrometric and RV data points. 

With a minimum mass of Mp sini = 2.98 Mj up , a period 
of P = 3.63 yr and an eccentricity ofe = 0.48, the solution 
for the giant planet (Orbit Ab) is consistent with the orbital 
parameters found by Correia et al. (2008) and Fischer et 
al. (2009). The imaging data lift the degeneracy of 
orbital parameters for the stellar companion HD 196885 B. 
The mass of 0.45 Mq is in agreement with the companion 
photometry and the M1±1V spectral type (Chauvin et al. 
2007). With a period of P = 72.1 yrs, our observations 
confirm a close orbit with a semi-major axis of just 21 AU. 
Together with G186 (Queloz et al. 2000; Lagrange et al. 
2006), 7 Cep (Hatzes et al. 2003; Neuhauser et al. 2007) 
and HD 41004 (Zucker et al. 2004), this system is among 
the closest binaries with one component hosting a giant 
planet. 

5. Dynamical evolution and stability 

We performed a preliminary numerical study of this system 
using the symplectic N-body package HJS (Beust 2003) op- 
timized for hierarchical multiple systems. The study was 
started from the fitted solution of Table [2] considering that 



Table 2. Fit and orbital parameters of HD196885 Ab and 
B 



Par am. 


[unit] 


HD 196885 Ab 


HD196885B 


P 


[yr] 


3.63 ±0.01 


72.06 ±4.59 


e 


0.48 ±0.02 


0.42 ±0.03 


0J 


[deg] 


93.2 ±3.0 


-118.1 ±3.1 


Q 


[deg] 




79.8 ±0.1 


i 


[deg] 




116.8 ±0.7 


t P 




2002.85 ± 0.02 


1985.59 ±0.39 


a 


[AU] 


2.6 ± 0.1 


21.00 ±0.86 


Maksim 


[M Jup ] 


2.98 ±0.05 




M B 


[Mq] 




0.45 ±0.01 


vo (Coravel) 


[km/s] 


-31.8S 


)± 0.20 


v (Elodie) 


[km/sj 


-31.99 ±0.04 


wo (Lick) 


[km/s] 


-1.64 


±0.04 


vo (Coralie) 

x 2 


[km/s] 


-32.00 ±0.04 




465.92 



some free parameters in the Ab orbit remained (namely the 
inclination i and the longitude of node £1). Depending on 
the inclination assumed, the mass of the giant planet could 
range between the minimum value quoted in Table [2] and 
several tens of Jupiter masses. The angle il had no effect 
on the RV signal of the planet, but is related to the mutual 
inclination i r (mutual inclination between the two orbits) 
by: 

cos i r = cos i cos i! ± sini sini' cos(S! — fl') , (1) 

i' and f2' are the corresponding parameters for the AB or- 
bit (listed in Table [2]) . i and i r turn out to be the more 
relevant parameters to dynamically characterize the sys- 
tem. We thus performed various integrations over 10 7 yr 
with different values of (i, i r ) couples. Not all pairs were 
compatible with our constraints. 

Two main results were obtained from our simulations. 
The first result shows that in all cases the system appears 
chaotic with significant (more or less erratic) changes of 
semi-major axis for the Ab orbit. The stability of the sys- 
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Fig. 4. Radial velocity observations and orbital solution of HD 196885 AB based on the CORAVEL, ELODIE, CORALIE 
and Lick surveys. Top Left, Radial velocity variation as a function of time with the best orbital solution overplotted 
presenting the complete set of RV observations. Top Right, Zoom-in on the RV variation due to the giant planet 
HD 196885 Aa superimposed to the RV drift due to the companion HD 196885 B. Bottom Left, Radial velocity variation 
of HD 196885 A due to the B component only as a function of the orbital phase. Bottom Right, Radial velocity variation 
of HD 196885 A due to the Ab planet only as a function of the orbital phase. 



tem needs to be therefore investigated over a much longer 
timescale. The second result shows that the mass of the gi- 
ant planet HD 196885 Ab (hence the Ab orbit inclination 
i) does not significantly affect the stability of the system. 
Additional simulations with similar i r but different i show 
similar behaviours. This can be explained by the dominant 
influence of the higher mass companion HD 196885 B as the 
giant planet behaves almost like a massless test particle in 
this system. The mutual inclination i r is however a key pa- 
rameter for the system stability. The behaviour of two dif- 
ferent solutions (low and high i r ) concerning the Ab orbit 
is shown on Fig. [5j The AB orbit remains slightly affected 
and is not considered here. In both cases, we note signifi- 
cant semi-major axis oscillations revealing a chaotic regime. 
The eccentricity oscillations are high, especially in the high 
i r case. In fact, all high i r configurations fall in the Kozai 
resonance regime (Kozai 1962). Under the effect of secular 
perturbations by the outer body, the giant planet is subject 
to a periodic evolution that drives it to a lower inclination 
but a very high eccentricity. This mechanism is active in 



non-coplanar hierarchical triple systems (Harrington 1968; 
Ford et al. 2000) and is therefore not surprising in our case. 
It was actually invoked to explain the high eccentricity of 
some extrasolar planets in binary systems (Mazeh 1997; 
Holman et al. 1997; Libert & Tsiganis 2009; Fabricky et al. 
2007). 

Surprisingly, our results show that the system is more 
stable in the Kozai regime. The Kozai resonance is an angu- 
lar momentum exchange process. In a pure 3-body system, 
it does not affect the semi-major axis evolution and the 
system stability. Therefore the real level of instability can 
be read from the evolution of the semi-major axis, which 
is more stable in the Kozai configuration (see Fig. [5]). This 
is a general trend present in all our simulations. For copla- 
nar systems, only marginally stable solutions are found. 
Consequently, our results suggest that a non-coplanar con- 
figuration, characterized by a Kozai regime, is more prob- 
able for HD 196885. The high eccentricity of the Ab orbit 
fit is another strong indication in favour of a Kozai regime. 
If confirmed, this system would constitute one of the most 
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Fig. 5. Examples of integration runs for the HD 196885 3 body system. The evolution of the inner orbit Ab is shown 
over 10 7 yr. Top, a typical coplanar solution, with semi-major axis (left) and eccentricity (right) evolution; Bottom, 
same for a typical highly i r solution with kozai oscillations (see eccentricity). 



compact non-coplanar systems known so far. Deeper inves- 
tigation will be needed to further constrain this issue and 
will be presented in a forthcoming paper. 

6. Discussion 

Stellar companions do influence the occurrence and the 
properties of planets in circumstellar orbit (i.e orbiting 
one of the binary compenent). Observations show that 
planets seem less frequent in close (< 100 AU) binaries 
(Eggenberger et al. 2007, 2008) and that massive short- 
period planets appear to be preferentially found orbiting 
one component of a multiple system (Zucker & Mazeh 2002; 
Eggenberger et al. 2004; Desidera & Barbieri 2007). Among 
the few tens of exoplanetary hosts, G186, 7 Cep, HD 41004 
and HD 196885 are of particular interest as they are short 
period binaries with semi-major axis lower than ~ 25 AU 
(see system properties in Table 3). In such close systems, 
we can wonder how their planets might have formed and 
survived to the close interaction with the outer binary com- 
panion. The presence of a binary companion is expected to 
truncate the protoplanetary disk, affect its eccentricity and 
impact the velocity dispersion and evolution of planetesi- 



mals which could lead to a planet formation hostile environ- 
ment (e.g. Thebault et al. 2006; Paardekooper et al. 2008; 
Xie et al. 2008; Marzari et al. 2009; Cieza et al. 2009). In the 
context of 7 Cep, Thebault et al. (2004) and Paardekooper 
et al. (2008) showed that core-accretion formation of a gi- 
ant planet was feasible, although it was probably slowed 
down (assuming coplanarity of the disk, planet and binary). 
Jang-Condell et al. (2008) found that disk instability was 
in addition more unlikely as a massive circumstellar disk 
and/or extremely high accretion rate as seen in FU Orionis 
events are required. Kley & Nelson (2008) then examined 
the dynamical evolution when a protoplanetary core has 
formed and begins the accretion phase. They constrained 
to (a p < 2.7 AU) the initial semi-major axis of the giant 
planet that could lead to a stable solution surviving inward 
migration and eccentricity variation. HD 196885 closely re- 
sembles the 7 Cep system. However, the higher mass ratio 
between the binary companion and the primary and the 
further location of the inner planet makes it even more 
challenging for core-accretion planet formation theories. 

Although coplanarity seems a reasonable assumption for 
the planet formation in close binaries (Hale 1994), our pre- 
liminary dynamical study shows that a non-coplanar config- 
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Table 3. Orbital parameters of small separation (a ~ 20 AU) binaries with one exoplanct around one component. 



Name 


Age 


Mi 


M B 


a B 


es 


M p 


a p 


ep 




(Gyr) 


(Mo) 


(Mo) 


(AU) 




(M Jup \ 


) (AU) 




HD 196885 


1.5-3.5 


1.3 


0.45 


21.0 


0.42 


3.0 


2.6 


0.48 


7 Cep 


6.6 


1.6 


0.41 


19.0 


0.41 


1.6 


2.0 


0.12 


G186 


2.0 


0.8 


0.45 


18.0" 


0.40 


4.0 


0.11 


0.05 


HD 41004 


1.6 


0.7 


0.40 + 0.02 6 


~ 20 b 


0.40 


2.5 


1.64 


0.39 



a predicted intial semi-major axis of 13 AU for the binary (Lagrange et al. 2006). 

b HD 41004 B is actually a SB1 binary composed of a M4V star and a Msb = 19 Mj up brown dwarf (Santos et al. 2002; Zucker 
et al. 2004). 



uration, characterized by a Kozai regime, is more probable 
for HD 196885. Its origin might therefore be questioned, as 
it seems difficult to build a non-coplanar system from a 
single disk. The giant planet has to form and survive in a 
non-coplanar system under the combination of the secular 
pertubation of the close binary companion and disk evolu- 
tion. This has to be tested. An alternative scenario could 
be that the binary was less compact in the past. Close 
encounters with other young stars could have led its or- 
bit to shrink (Malmberg et al. 2008). The companion may 
also have been captured after the formation of giant planet 
around the primary in the birth cluster (Portegies Zwart 
& McMillan 2005; Pfahl & Muterspaugh 2006). We could 
even think of a capture of the giant planet itself. These are 
of course speculations but the high present-day eccentricity 
of HD 196885 B is also an indication of a chaotic past dy- 
namical history of that system. Further investigation will 
be needed to further constrain and test the system stability 
and origin on a broader parameters range. 

7. Conclusion 

We have reported the results of four years of astrometric 
monitoring of the planet-host binary system HD 196885 AB 
using NaCo at VLT. Combined with RV observations, our 
imaging results enabled to derive the inclination, the eccen- 
tricity and the true mass of the B component, a low mass 
star orbiting at 21 AU. We found consistent solutions for 
the planet with results previously reported in the litera- 
ture. We also confirmed that HD 196885 AB belongs to the 
rare cases of close binaries with one component hosting a 
giant planet, offering an ideal labolatory to study the for- 
mation and evolution processes in such extreme planetary 
systems. Finally, we have run N-body numerical simulations 
to test the system stability. The main result suggests that 
non-coplanarity and high mutual inclination, characterized 
by a Kozai regime, favor the system stability. If confirmed, 
HD 196885 would constitute one of the most compact non- 
coplanar system hosting a circumstellar planetary system. 
How planet formation could have occured in such an ex- 
treme and hostile environment remains a challenging ques- 
tion. In-situ formation of a non-coplanar system cannot be 
excluded and has to be tested. However, alternative sce- 
narii such as an external perturbation with stellar encoun- 
ters modifying the binary properties or capture mechanisms 
are likely to offer reasonable explanations. 
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